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ABSTRACT 

Context. The mass accretion rate (M acc ) is a key parameter that has not accurately been determined for a wide sample of Herbig Ae/Be 
(HAeBe) stars until recently. 

Aims. We look for trends relating M acc and the stellar ages (f), spectral energy distributions (SEDs), and disk masses for a sample of 
38 HAeBe stars, comparing them to analogous correlations found for classical T Tauri stars. Our goal is ultimately to shed light on 
the timescale and physical processes that drive evolution of intermediate-mass pre-main sequence objects. 

Methods. Mass accretion rates obtained by us in a previous work were related to several stellar and disk parameters: the age of the 
stars was compiled from the literature, SEDs were classified according to their shape and the wavelength where the IR excess starts, 
near- and mid-IR colour excesses were computed, and disk masses were estimated from mm fluxes. 

Results. M acc decreases with the stellar age, showing a dissipation timescale r = 1.3*1'? Myr from an exponential law fit, while a 
power law yields M acc (t) oc r' ; , with n = 1.8*i'i. This result is based on our whole HAeBe sample (1-6 M Q ), but the accretion rate 
decline most probably depends on smaller stellar mass bins. The near-IR excess is higher and starts at shorter wavelengths (J and H 
bands) for the strongest accretors. Active and passive disks are roughly divided by ~ 2 x 10~ 7 Mq yr _1 . The mid-IR excess and the 
SED shape from the Meeus et al. classification are not correlated with M acc . Concerning disk masses, we find M acc oc M disk u±0 - 3 . Most 
stars in our sample with signs of inner dust dissipation typically show accretion rates ten times lower and disk masses three times 
smaller than the remaining objects. 

Conclusions. The trends relating M acc with the near-IR excess and M djsk extend those found for T Tauri stars, and are consistent 
with viscous disk models. The differences in the inner gas dissipation timescale, and the relative position of the stars with signs of 
inner dust clearing in the M acc -M rfl0 >. plane, could be suggesting a slightly faster evolution, and that a different process - such as 
photoevaporation - plays a more relevant role in dissipating disks in the HAeBe regime compared to T Tauri stars. Our conclusions 
must consider the mismatch between the disk mass estimates from mm fluxes and the disk mass estimates from accretion, which we 
also find in HAeBe stars. 

Key words. Stars: pre-main sequence - Stars: variables: T Tauri, Herbig Ae/Be - Accretion, accretion disks - circumstellar matter - 
protoplanetary disks 

1 . Introduction could in turn be estimated from the amount of gas accr eted by 

the star (lAlexander & Armitagell2007tlNaiita et al.ll2007l) . 

Nowadays there is cons ensus that practica lly all stars are formed 

from accretion disks dKraus et all |20TqT) , even if these have Mass accretion rates in classical T Tauri (CTT) stars are 
sometimes very short dissipation timescales. The mass accre- mainly derived from the UV emission excess and spectroscopic 
tion rate (M acc ) is a major parameter that drives the star-disk veilin g caus ed by the accretion shocks, as well as from mag- 
evolution, given that a significant fraction of the stellar mass ac- netospheric accretion line modelling and empirical correlations 
quired during the pre-main sequence (PMS) phase is accreted with the luminosity of several emission lines (see the review by 
from the disk. Circumstellar disks are also dissipated through ICalvet et al|2000j and references therein). The CTTs show a de- 
other physical processes, including photo-evaporation, dust set- creas e in M acc of about one or two orders of magnitude in the first 
fling and dynamical interactions with forming planets (see e.g. M Y r (e.g.|Fang et aTJ|2009|), which gives an estimate of the inner 
ICiezal l2008t I Williams & Ciezal |20T1 . Their relative influence g as dissipation ti mescale. This is intim ately related to the inner 

dust disk lifetime (Hart igan et alJ[l995l) . indicating that both gas 

and dust close to the star vanish at a similar rate. Measurements 

Send offprint requests to: Ignacio Mendigutia of dust emission at longer wavelengths provide estimates of the 

e-mail: Ignacio.Mendigutia@cab.inta-csic.es disk mass (M^) and point to an inside-out dissipation of the 
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circumstellar environment (Cieza 2008, and references therein). 
The so-called transitional disks, where the IR emission is weak 
or negligible up to A > 10 /mi, are particularly interesting. Stars 
with this type of spectral energy distribution (SED) are inter- 
preted to be in a fast evolutionary stage that bridges the gap be- 
tween the optically thick disk phase and the optic ally thin debris 
disk one (Wol k & Walterlll99rl ICieza et al.ll2007l) . An analysis 
of the Mdisk-M aC c plane can be useful to discern between the 
different physical processes that drive disk dissipation. In partic- 
ular, the region occupied by several CTTs with transitional disks 
in that plane is com patible with the formation of Jovian planets 
(Naiit aetal.ll2007l) . 

Herbig Ae/Be (HAeBe) objects dWatersi |2006) are the 
intermediate-mass (~ 1-10 Mq) counterparts of CTT stars. The 
HAeBes are the most massive objects to experience an opti- 
cally visible pre-main sequence phase, bridging the transition 
between low- and high-mass stars. Although these are highly in- 
teresting objects - e.g., planet formatio n could be more efficient 
in stars more massive than the Sun dKennedv & Kenvonll2008t 
Bow ler et al.ll2010i: lBossll20lTh - our knowledge about the prop- 
erties of the HAeBe regime is much more limited than that for 
the T Tauri stars. The main reason is the comparatively smaller 
HAeBe sample, caused by the faster evolution of massive stars 
to the main-sequence, and by the fact that star formation favours 
lower mass objects, as the shape of the initial-mass function sug- 
gests (lSalpeterlll995l) . 

The optical photometry in Oudmai ier et al.l (1200 ll) and the 
simultaneous optical spectroscopy in Mendigutfa et ajj (1201 laL 
Paper I hereafter) were recently used to derive accretion rate esti- 
mates and_ejn2irj£alrelati£ns^ with emission lines for 38 HAeBe 
stars (Mendig utfa et al.l 1201 lbi Paper II hereafter). The accre- 
tion rates in Paper II are the first obtained for a wide sample of 
HAeBe stars from a specific modelling of the UV B aimer ex- 
cess and constitute the most reliable estimates for this regime 
to date. The present work closes this series of papers. Here, we 
analyse possible relationships between the mass accretion rates 
and several stellar and disk parameters of the 38 HAeBes stud- 
ied in Papers I and II, and compare them to the better known 
properties for CTTs. In particular, we look for trends with the 
stellar age, properties of the SED, and disk mass. The ultimate 
objective is to contribute to the knowledge of circumstellar disk 
dissipation in HAeBe stars. The paper is organized as follows. 
Section [2] describes some general properties of the sample. The 
SED properties and disk masses are treated separately in Sects. 
I2.1l and [272l respectively. The correlations with the mass accre- 
tion rate are presented and analysed in Sect. [3] The discussion is 
included in Sect. [4] and the summary and conclusions in Sect. [5] 

2. Sample properties and data 

The sample consists of 38 HAeBe stars covering almost all 
such obje c ts in t he northern hemisphere from the catalogue of 
iThe et al.l (11994 (see Papers I and II). All stars show IR ex- 
cess and variable emission lines (Merfn 2004, Paper 1). Columns 
2 to 8 in Table [1] include the effective temperatures, distances, 
ages, mass accretion rates, and information on possible stellar 
companions. The parameters in Cols. 9 to 16 are described in 
the following sections. The spectral types cover the Ae and late- 
type Be reg ime, including ten intermediate-mass F and G stars 
dMoraet al.ll2001l) . The stellar age covers the pre-main sequence 
phase for our sample, up to ~ 15 Myr. Most ages were deter- 
mined by Mon tesinos et alJ (|2009), with a typical uncertainty of 
35%. The range in stellar mass and mass accretion rate is ~ 1-6 
M and 10 9 — 10 5 M Q yr , with a typical -median- value of 2.4 



X 10 1 M yr -1 . The mass accretion rates are those reproduc- 
ing the observed Balmer excesses from a magnetospheric accre- 
tion shock model (Paper II). The accretion rates for R Mon, VY 
Mon, VV Ser, and LkHa 234 were estimated extrapolating the 
calibrations with the emission line luminosities provided in that 
work. Their large uncertainties reflect the difficulty to model the 
strong Balmer excess of these stars. In case of multiplicity, most 
of the stellar and accretion parameters are in principle referred 
to the brightest component, but some contamination cannot be 
excluded for objects with close companions (d < 1"; see refer- 
ences in Table [TJ. T he high percentage of binaries is typical for 
HAeBe stars dWheelwright et al.ll2010b . 

2.1. Spectral energy distributions 

Multi-wavelength photometry compiled from the literature is 
listed in Tables |A~T1 IAT21 and |A~31 Figure |A~T1 shows the SEDs 
of the stars in the sample. The best fits obtained using the on- 
line SED fitting tool f r om th e 2-D radiative transfer models in 
iRobitaille etaiT (I2006L I2007IF1 are overplotted. The use of si- 
multaneous U BVRIJ HK photometry se l ected a t the brigh test V 
state f rom the data in lOudmaiier et al.l (12001 1) and lEiroa et all 
(2002) allow us to minimize the influence of obscuring ef- 
fects from circumstellar material, providing the best possible 
photospheric-innerdisk fitting. The wavelength where the IR ex- 
cess becomes apparent can thus be derived for each object (Sect. 
12.1.21 ). Photometry at wavelengths longer than 160 //m was not 
used if its inclusion made the stellar and inner-mid disk fitting 
worse - note that Robitaille's models were originally optimized 
for the 1 to 100 /urn region -. The stars are classified according to 
the shape of their SEDs and the properties of their IR excesses, 
as it is described in the two following sections. 

2.1.1. Classification from the SED shape 

Colu mn 9 in Ta ble Q] shows the SED group according to 
the iMeeus et al.l ((2001) classification scheme. The Meeus 
group I and group II (M01 groups hereafter) SEDs have 
been associated with geometrical differences in the struc- 
ture of the disks , flared and self-shadowed, respectively 
dMeeus et alJ 12001b iDullemondl [20021: iDullemond & Dominikl 
2004). Both groups s how differences regardin g the photopo- 
larimetric behaviour (IDullemond et al.l |2003), presence and 
streng th of the polycyclic aromatic hydrocarbons (PAHs) fea- 
tures (lAcke & van den Anckei 20041). and gra in growth as ob- 
served from sub-mm data ( Acke et ail 120041) . The classifica- 
tio n in M01 groups was m ade following the colour criterium 
in Ivan Boekel et al.l d2003l) . This is based on the ratio L(nIR) 
(the integrated luminosity obtained from broad-band JHKLM 
photometry) to L(IR) (the equivalent parameter, obtained from 
the IRAS 12, 25 and 60 micron fluxes), and on the non- 
colour corrected IRAS [12] - [60] colour. This classification 
used the photometry given in Tables I A. 1 1 and IA.2I and agrees 
with previous one s reported for many stars in our sample 
(IMeeus et all 120011: lAcke et all 120041; lAcke & van den Anckerl 
12004 1 Acke et alJl2010h . A direct visual inspection of Fig. |AJ] 
allows one to classify HD 58647, BH Cep and NV Ori as M01 
group II sources, and AS 442, and VY Mon as group I. T Ori is 
provisionally classified as a group I source, but mid-IR photom- 
etry is necessary to confirm this result. Summarizing, 12 sources 
in our sample are classified as group I stars, the remaining 26 
being group II objects. 



http://caravan.astro.wisc.edu/protostars/ 
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Star 


T 


d 


Age 


Mace 


d(SpT') 


SK 


Ref 


M01 




E(H - K) 


[12]-[25] [25]-[60] 


T 


P 


M disk 




(K) 


(pc) 


(Myr) 


CM© yr 1 ] 


(") 


(mag) 






(band) 


(mag) 


(mag) 


(mag) 


(K) 




[Mo] 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


HD 31648 


8250 


146 


6.7 


<-7.23 








II 


H 


0.71 


-0.18 


0.08 


21 


0.9 


-1.27+0.03 


HD 34282 


9550 4 


164" 


>7.8" 


<-8.30 








I 


K 


0.57 


0.92 


2.01 


38 


1.4 


-1.94+0.04 


HD 34700 


6000 s 


336 £ 


3.4 s 


<-8.30 


5.2(Ml-2) 


3.3 


(1) 


I 


>K 


0.18 


2.17 


1.26 


16 


1.5 


-1.86+0.04 












9.2(M2-3) 


4.2 


(1) 


















HD 58647 


10500 


543 


0.4 


-4.84±0.22 


>0.5(?) 


9 


(2) 


11 


H 


0.65 


-0.59 


-1.96 








HD 141569 


9550 4 


99 A 


6.7" 


-6.89±0.40 


7.6(M2) 


1.8" 


(3) 


II 


>K 


0.04 


1.33 


1.18 


33 


1.2 


-3.41+0.19 












9.0(M4) 


2.4* 1 


(4) 


















HD 142666 


7590 4 


145 A 


5.1" 


-6.73±0.26 








11 


K 


0.56 


0.29 


-0.48 


21 


0.6 


-1.73+0.03 


HD 144432 


7410 4 


145" 


5.3" 


<-7.22 


1.4(K5Ve) 


2.4 


(4,5) 


11 


H 


0.61 


0.24 


-0.53 


16 


1.0 


-2.15+0.04 


HD 150193 


8970 


203 


5.0 


-6.12±0.14 


l.l(F9Ve) 


2.2 


(4,6) 


11 


J 


0.67 


0.03 


-0.87 


26 


0.2 


-1.98+0.04 


HD 163296 


9250 


130 


5.0 


<-7.52 








II 


H 


0.76 


0.15 


0.32 


29 


0.5 


-1.18+0.02 


HD 179218 


9500 


201 


3.3 


<-7.30 


2.5(?) 


6.6 


(7,8) 


I 


K 


0.64 


0.67 


-0.41 


29 


1.3 


-1.84+0.04 


HD 190073 


9500 


767 


0.6 


-5.00+0.25 








II 


H 


0.74 


-0.28 


-1.15 








AS 442 


11000 


826 


1.5 


-5.08±0.11 


?(?) 


? 


(9) 


I 


I 


0.86 


-0.71 


1.38 








VXCas 


10000 


619 


6.4 


-6.44±0.22 


5.3(?) 


4.8 


(8) 


II 


K 


0.91 


0.59 


-0.69 


28 


1.0 


<-1.92 


BHCep 


6460 4 


450* 


8.2" 


<-8.30 








11 


K 


0.55 


-1.53 


0.19 








BOCep 


6610 4 


400* 


11.2" 


-6.93+0.28 








I 


> K 


0.19 


1.52 


1.75 








SVCep 


10250 


596 


5.2 


-6.30+0.20 








II 


H 


0.72 


0.23 


-0.73 


24 


1.0 


-1.82+0.12 


VI 686 Cyg 


6170 4 


980" 


<0.2" 


-5.23+0.41 








I 


K 


1.05 


0.53 


1.50 


16 


1.0 


<-0.42 


RMon 


12020" 


800" 


< 0.01" 


(-4+1) 


0.7(K1V) 


6> 2 


(10) 


I 


H 


1.49 


0.96 


-0.09 


56 


0.7 


-0.78+0.04 


VY Mon 


12020" 


800" 


< 0.01" 


(-4+1) 








I 


I 


0.97 


0.67 


0.57 


32 


2.4 


-0.45+0.04 


51 Oph 


10250 


142 


0.7 


<-7.00 








11 


>K 


0.42 


-0.47 


-2.46 








KK Oph 


7590 4 


160" 


3.9" 


-5.85+0.15 


1.5(G6Ve) 


2.5 


(4,6) 


II 


K 


1.33 


-0.03 


-0.48 


16 


0.1 


-1.91+0.04 


TOri 


9750 


472 


4.0 


-6.58+0.40 


7.7 


>4.5 


(11) 


I 


>K 


0.86 






30 


1.0 


-0.85+0.09 


BFOri 


8970 


603 


3.2 


<-8.00 








II 


K 


0.60 


-0.20 


1.22 


24 


1.0 


-1.88+0.14 


COOri 


63 10" 


450" 


< 0.1" 


-5.20+0.18 


2.0(?) 


3.2* 3 


(12) 


II 


H 


0.65 


-0.13 


0.40 


15 


1.0 


<-2.11 


HK Ori 


85 10 4 


460" 


1.0" 


-5.24+0.12 


0.4(K3) 


2.30 


(11,13) 


II 


>K 


0.86 


0.08 


-0.99 


20 


1.0 


<-0.81 


NVOri 


6750 c 


450 F 


4.4 C 


<-8.30 








II 


K 


0.51 












RYOri 


6310" 


460" 


1.8" 


-6.65+0.33 








II 


H 


0.49 


-0.06 


1.27 








UX Ori 


8460 


517 


4.5 


<-6.77 


> 0.02(?) 


7 


(12) 


II 


H 


0.72 


0.35 


-0.28 


28 


0.6 


-1.49+0.04 


V346 Ori 


9750 


586 


3.5 


-5.99+0.17 








I 


K 


0.47 


1.51 


1.63 








V350 Ori 


8970 


735 


5.5 


-6.66+0.24 


0.3(?) 


3.2 


(8) 


11 


>K 


0.71 


0.67 


-0.42 








XY Per 


9750 


347 


2.5 


-5.86+0.20 


1.2(7) 


0.0 


(6) 


II 


H 


0.73 


0.06 


0.21 








VVSer 


13800 


614 


1.2 


(-5.2+0.8) 








II 


H 


0.93 


-0.22 


0.56 


45 


1.0 


<-1.61 


CQTau 


6800 s 


130° 


7.7 s 


<-8.30 








I 


K 


0.73 


1.31 


0.06 


14 


0.6 


-1.60+0.03 


RRTau 


10000 


2103 


0.4 


-4.11+0.16 


?(?) 


? 


(7) 


11 


H 


0.88 


0.25 


0.77 


28 


1.0 


<-0.33 


RYTau 


5770 c 


134 F 


6.5° 


-7.65+0.17 








II 


H 


0.68 


0.43 


-0.58 


14 


1.0 


-1.50+0.04 


PX Vul 


6760 4 


420" 


14" 


-6.72+0.16 








II 


K 


0.57 


0.49 


-0.68 








WW Vul 


8970 


696 


3.7 


-6.38+0.70 








11 


H 


0.83 


0.31 


-0.38 


24 


1.0 


-1.51+0.04 


LkHa 234 


12900" 


1250" 


< 0.01" 


(-4.1+0.4) 


2.7(7) 


>5.20 (11) 


I 


H 


0.88 


1.82 


2.35 


65 


2.5 


0.31+0.04 



Notes. Columns 1 to 5 list the name of the star, stellar effective temperatu re, dis tance, age, and mass accretion rate. Unless otherwise indicated, the 
stellar temperatures, distances and ages are taken from Montesinos et al. (2009). Accretion rates (on log scale) are taken from Paper II. Columns 6 
to 8 refer to possible stellar companions. Column 6 lists their angular distance and spectral type, Col. 7 their difference in K magnitude (K s for *', 
K' for * 2 and Hipparcos H ; , for * 3 ), and Col. 8 the corresponding references. HD 34700 is a triple visual system, the brightest object of which is a 
spectroscopic binary (1). HD 141569 is a triple system (3). Columns 9 to 13 list the SED classification from the Me eus et al] d200lh scheme, the 
photometric band according to the shortest wavelength where the IR excess is apparent, the intrinsic near-IR colour excess, and the IRAS colours. 
Columns 14 to 16 list the dust temperature, /? value, and disk mass (on log scale). These were obtained from the fluxes at 1.3 mm, except for HD 
34700, V168 6 Cyg, RY Ta u and L k Ha 234 (see TablelA2ll. 

References. iManoi etaD j2006h. 'lAlonso-Albi et al] d20"09|) and refer e nces t herein, HMerinl d2004h. ^Siess et al] dl999h ^Acke etafl d2005l) . 
iBlondel et al] "" 
Carmonaj 



HMerinl d2004h. ^ISiess et al.1 dl999h ^Acke etal.1 

. _ _ . d2006h . ( 1): ISterzik et al] d2005l) . ( 2): iBaines et al] d2006D, ( 3): IWeinberger et al] EoOO). (4): 

Perez et all d2004h. (6) : iPirzkal et al.l dl997h. (7 ): IWheelwright et all d20ldh . ( 8): iThomas et all d2007l) . (9): 



iGarcia Lopez et al 
(5): 



Corp oron & Lag range ( 1999b. (10): IClose et al.l(il997l) , (ll): lLeinert et al.ldl997l) . (12): iBertout et al] d 19991) . (1 3): ISmith et al] d2005) 



2.1.2. IR excess 

The sample was divided into two groups according to the short- 
est wavelength where the IR excess is apparent from the SED 
fits in Fig. IA.1I Half of the stars show an IR excess starting at 
wave lengths correspondi ng to the J or H bands (1.22 and 1.63 
//m; iBesse ll & Brettl l 1988h . and are included in the first group 
(group JH hereafter). The objects in the second group (group 
K) show an IR excess starting at A > 2.19 //m (filter K). The dust 
in the disks of group JH sources is expected to extend up to the 



sublimation radius (around 0.5 AU, for a typical HAe star with 
T„ = 9000 K and R* = 3 R fl Assuming that the dis k temper- 
ature radially declines following a R~ 3 ^ 4 law (see e.g. lArm itage 
2009, and references therein), the small dust grains in group K 
sources would be located at least twice far away from the central 
star. In particular, the SED of the four stars that show transition 



2 The dust sublimation radius is R s „j, ~ R»(Tj„;,/T„)~ 2 , with R„, T« 
and T IBi , the stellar radius, effecti ve temperature and s ublimation tem- 
perature (~ 1600 K), respectively dRobitaill e et al. 2003). 
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disks - defined as showing very weak or null excess up to ~ 10 
fan: HD 34700, HD 141569, BO Cep, and 51 Oph - could be the 
signature of inner holes with sizes larger than the dust sublima- 
tion radius by an order of magnitude. Given that R SM j increases 
with the temperature and radius of the central object, it is worth 
noting that the possible inner holes associated to HAeBe stars 
with IR excesses starting at the K band could in principle have 
sizes similar to tho se associate d to transitional T Tauri stars (~ 
a few AU; see e.g.|Merin 2010). Alternatively, the region in the 
disk sampled by an IR excess starting at a given wavelength is 
farther away from the central source in HAeBe stars than in TTs. 
The photometric band where the IR excess starts for each object 
is shown in Col. 10 of TableQ] 

Almost all stars in the sample have measurements in H, K 
and the IRAS bands, whose emission will be used to probe the 
disk at different radial distances. Column 11 in Table Q] shows 
the intrinsic colour excess E(H - K)q = (H - K)d er ed - (H - K)q. 
This parameter m easures the colour excess without the extinc- 
tion contribution dMever et al.l 1 19971) . reflecting the properties 
of the inner disk dust. The de-redenned colour, (H - K)dered 
= (H - K) + R V E(B - V) x (k/r/ky - k H /ky), k^ being the 
opacity at a given wavelength, is derived considering the extinc- 
tion law in Robitail fe et al.l (120071). a to tal-to-selective extinction 
ratio Ry = 5 dHernandez et alj 120041) . and the B-V colour ex- 
cesses from Paper II . The in trinsic colours (H - K)o are from 
iKenvon & Hartmannl dl995h . Columns 12 and 13 of Table [Q 
show the IRAS [12]-[25] and [25]-[60] colours, reflecting the 
warm mid-disk dust. The uncertainty in the IRAS colours - typ- 
ically 0.15 magnitudes - is significantly larg er than the reden- 
ning and intrinsic colour corrections (see e.g. iPatten & Willsonl 
Il99ll) . which were therefore not applied for these bands. 



2.2. Disk masses 

Column 16 in Table Q] shows the total (dust + gas) disk masses 
for the stars with mm photometry (see Table IA.3}. Estimates o f 
the disk masses were derived following Beckwith et al ] (fl990l) . 
assuming that the emission is optically thin at mm wave- 
lengths. The following expression, valid under the Rayleigh- 
Jeans regime, has been used: 



have one single mm measurement. For these objects we followed 
Natta et all (1200 0). fixing B to 1 and To to the corresponding 
value from the spectral type, according to their Table II. The fi- 
nal values for To and B are listed in Cols. 14 and 15 of Table 

M disk ranges typically between 10 4 and 0.4 Mq, 0.02 Mq 
being the median value. Considering the stellar masses, the me- 
dian value for the M^/M* ratio is 1 % (i.e. the same as for 
lower-mass T Tauri stars; see lAndrews & W illiams 2007). The 
extremely high disk mass for LkHa 234 comes from the flux 
data set in Pezzuto & Strafella] (1 19971) (see Table IA. 31 . which is 



probably contaminated by environment emission. iFuente et al.l 
(2001) reported a flux of less than 20 mJy at 1.3 mm using a 
smaller 1.2" x 1.3" beam, deriving a disk mass < 0.1 M for 
LkHa 234, whereas lAlonso-Albi et alj d2009l) reported M*-,* < 
0.2 M from the same flux. It is finally noted that the uncertain- 
ties for Mdisk listed in Table Q] only consider the errors in the 
mm fluxes. Therefore, they should be considered lower limits, 
since additional error sources are most probably coming from 
the Td and B values used in the fits, and specially from the un- 
certainty in the opacity (see Section |4] and references therein). 
Gas-t o-dust ratios significantly lower t han the standard va lue of 
100 (iPanic et alJl2008HThi et alJl2010l:lTilling etai1l2012h could 
also affect the disk masses derived. For instance, our disk masses 
are an order of magnitude larger than those obtained from Eq. Q] 
and a gas-to-dust ratio of 10. 



3. Correlations analysis 

In this section we analyse possible correlations between the mass 
accretion rates from Paper II (Col. 5 in Table[T) and several prop- 
erties of the stars in the sample, described in previous sections. 



M, 



disk 



2kT D v 2 k v 



(1) 



where F v is the flux measured at 1.3 mm for most stars in the 
sample, d is the distance, k is the Boltzmann constant, Td is 
the dust temperature, and k v the opacity at the observed fre- 
quency. We adopted the standard ki.3 mm = 0.02 cirr 2 gr -1 (as- 
suming a gas-to-du st ratio of lOOf j and a wavelength depen 



dence k mm oc A 13 (Beckwith et al. 



The value of B is representative of the maximum dust grain size, 
which is ~ 2-2.5 for the IS medium and sm aller for larger grains 
dPollack & Hollenbachl[T994l:lDrainell2006l) . The parameters T D 
and B are those best fitting t he available fluxes longward 35 
fim from a graybody (see e.g. lAndre et at.lll993l [Sandell 2000). 
The input range for the dust temperatures was within ~ ± 15 
K around a central value taken from the relation between this 
parameter and the spectral type (Natta et al.1 l2000l) . Figure IA.2I 
shows the results from these fits. The typical B value obtained 
from them is ~ 1.0 (mean and median), pointing to grain growth 
with respect the interstellar medium. Several stars in the sample 



3. 1 . Accretion rate and stellar age 

The mass accretion rates are plotted against the stellar ages in 
Fig-E The accretion rate declines as the age increases, following 
a trend that can be fitted to an exponential expression M. dcc {f) = 
Micc(0) x exp (—t/r), being t the stellar age and t the dissipation 
timescale (Ma noi et alJl2006tlFedele et al.ll2010t) . Apart from in- 
trinsic scatter, the final value for t is affected by uncertainties in 
the accretion rates and stellar ages, and also by upper and lower 
limits in both parameters for several stars. Both sources of er- 
ror can be treated statistically in independent ways, i.e., there 
is no statistical technique able to deal simultaneosly with un- 
certainties for detections and uppe r/lower limi t estim ates, to our 
knowledge. We used the method in [York et al .1 (120041) to derive a 
fit assuming different ranges of reasonable uncertainties for both 

1990) with < B <2 5 ^ acc anc ^ f ' obtaining T = 0-9^ 2 Myx. On the othe r hand, the 

astron omi cal survival ana lysis methods discussed in Isob e et ail 
( 1986) and lLavallev etafl l l 19921) are useful to deal with censored 
data that include upper and lower limits. The schmidttbin task 
in the iraf/ stsdas package is based on those works and pro- 
vides a value of t = 1.7 ± 0.6 Myr. Taking the average from 
both estimates, we assumed t = 1.3 
dissipation timescale for our sample 

The accretion rate decrease with the stellar age can also be 
fitted by a power-law express ion M acc (t) oc t~ n , wh ich is used 
in viscous dissipation models (Hartm ann et al.l fl998). From the 
statistical methods explained above, our data can be reproduced 



_o 5 Mv r as the accretion rate 



3 The total (dust + gas) opacity is k 13mra = kf^/fr+l), with k*^ im 



the dust opacity and r the disk mass gas-to-dust ratio. 



with rj = 2.3 | q, taking uncertainties for the ages and accretion 
rates into account. If the upper and lower limits are considered, rj 
— 1.2 + 0.1. The average fit for our sample is rj — 1.8*il which 
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Fig. 1. Mass accretion rates versus stellar age for our HAeBes 
(solid symbols, upper and low er limits are represente d with tri- 
angles) and T Tauri stars from lHartmann et al.l (11998b (crosses). 
The best fits for our sample are represented with solid and 
dashed lines, for the exponential (r = 1.3*i'? Myr) and power- 
law (77 = 1.8*H) expressions, respectively. 51 Oph was not con- 
sidered for the fit. The dotted line is the best power-law fit for T 
Tauri stars (77 -1.5; se e lHartmann et alJl998tlCalvet et alJ200Ct 
iMuzerolle et"aT1l2000l) . 

is shown in FigQ] along with the best fit obtained for CT T stars 
with K and M spectral types from lHartmann et al] d 1998b . 

The exponential decay provides a better fit than the power- 
law one in terms of a lower x 2 , and could reflect better the high 
but finite accretion rates of the youngest sources in the sample. 
However, this conclusion depends on the strong uncertainties for 
the accretion rates and ages of the youngest stars, which are also 
the most scarce in our sample. The stellar a ges for R Mon, VY 
Mon, and LkHa 234 are lower than 10 4 yr (iManoj et al.l l2006). 
and were derived from a different meth odology and evolutionar y 
tracks than for most star in the sample (Mo ntesinos et al.l l2009). 
Those values may be not accurate, but they reflect the youth- 
ness of these objects, as is supported by othe r estimates, which 
provide ages < 1 Myr for these stars (see e.g.[AJonso-Albi et al. 
2009; Tetzlaff et alJl20llT) . Our exponential fit provides a longer 
dissipation timescale than if their ages were older than assumed 
here. In this sense, the value for r obtained in this work could 
be considered as an upper limit. Analogously, the power-law fit 
would be steeper if the objects mentioned are older. In addition, 
5 1 Oph stands apart from the trend, and was not considered in 
the fits. The evolutionary stage of this object is not clear, show- 
ing characteristics typical of both a HB e star and a MS star with 
a gas -rich debris disk (Ivan den Ancker et al. I l200Tt IStarketaD 
2009, Paper II). 

The exponential expression that fits the accretion rate decline 
with age for our sample provides a timescale shorter than that ob- 
tained for the Ha equivalent width of HAeBe stars (~ 3 Myr in 
IManoi et al.l2006l) . However, the equivalent width itself does not 
measure the gas content and is subject to st rong variability (see 
Paper I). The difference with the estimate bv lManoi et al.l (|2006) 
can be diminished by transforming the line equivalent widths 



into line luminosities, and then translating these into accretion 
rates (See Paper I and II). The CTT sta rs show a decrea se of ~ 
1-2 dex in log M acc during the first Myr (Fan g et al.l2009l) . as our 
sample roughly show. Our value for t is lower than t he ~ 2.3 Myr 
report ed for K0-M5 pre-main sequence stars by iFedele et al.l 
(2010), although they are comparable considering the large un- 
certainties. Those authors made a statistical study of the objects 
that show signs of accretion in star forming regions at different 
ages, concluding that planet formation and/or migration in the 
inner disk might be a viable mechanism to halt further accre- 
tion onto the cent ral star on such a short timescale. Indeed, c ore- 
accretion models (Polla ck et al.l[l9 96; Mordasini et al. 200 ^) re- 
quire ~ 3 Myr including planet migration (Alibert et al. 2004, 
120051) . while gravitational instability dDurisenl2007l: lBosslf20lTE 
needs 10 3 -10 6 yr to form planets. The fast inner gas dissipation 
timescale suggested from our sample is more consistent with 
gravitational instability models, and contrasts with that required 
for the formation of planets from core-accretion ones. Other pro- 
cesses apart from planet formation might contribute to explain 
the fast decrease of the accretion rate in our stars, as is discussed 
in following sections. 

The analysis of the accretion rate decline represented as a 
power-law also points to a slightly faster inner disk gas dissi- 
pation in HAeBes, given that the mean 77 value derived for our 
sample is higher than that for CTT stars. Despite the strong 
unce rtainties, 77 ~ 1.5 seems more consistent for low-mass ob- 
jects dHartmann et alj|1998t ICalvet et alj|2000t IMuzerolle et~aTI 
2000), and the most recent estimates indicate an even slower 
mass accretion rate decay ( 77 ~ 1 .2: ISicilia-Aguilar et alJ bOlO: 
ICaratti o Garatti et al.ll2012b . 

We have followed a similar procedure as used for CTT stars 
by directly comparing the M acc and t values in the HAeBe 
regime. In this way, we have shown that the accretion rate and 
the stellar age are strongly correlated in our sample (Fig. [0, with 
a Spearman probability of false correlation of only 8.5 xlO 6 
%. However, the interpretation above of this trend also has to 
consider a caveat that comes from internal dependences on the 
stellar mass. First, M acc depends on it as ~ M, for our sam- 
ple (Paper II). Second, the age of HAeBe stars is also corre- 
lated with M t , the you ngest objects being the most massive 
(Ivan Boekel et alJl2005l Paper II). The partial correlations tech- 
nique (e.g. I Wa ll & Jenkins 2003) provides the probability of 
false correlation between two variables, with the effect of a third 
controlling variable removed. The Spearman probability of false 
correlation between M acc and t for our sample in Fig. [1] rises 
to 85 %, if the dependences of both M acc and t on M* are con- 
sidered. As a consequence, there is a high probability that the 
accretion rate decline in Fig. |T] is driven by the dependence on 
the stellar mass. Therefore, the r and 77 values describing this 
evolution should be viewed with caution. These results suggest 
that the evolution of the accretion rate should be studied by di- 
viding the HAeBe regime into small stellar mass bins, which 
would therefore need larger samples than studied here. This 
division would also be physically consistent, given that differ- 
ent evidence suggests that the accretion paradigm could change 
at some point within the HAeBe regime, the less massive ob- 
jects showing properties more similar to CTT stars (IVink et al.l 
120021 lEisner et al.ll2"004t IMottram et afl2007l Paper I, Paper II). 
Almost all stars in our sample with M, > 2.5 M are younger 
than 1 Myr, which could indicate that at older ages most of these 
objects have already reached the MS phase. The faster evolu- 
tion of the massive HAeB es would agree with recen t observa- 
tional works; in particular, iRoccatagliata et all (1201 lb reported 
that sources with masses > 2 M in the OB association IC 1795 
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Fig. 2. Histogram showing the mass accretion rates for the 38 
stars in the sample (shaded bars), and the stars belonging to 
the groups I and II (12 an d 26 objects ; solid and dashed lines) 
from the classification of iMeeus et al. | d200Tb . The probability 
that both groups are drawn from the same accretion rate distri- 
bution is 0.75. 
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Fig. 3. IRAS colours [12]-[25] (open symbols) and [25]-[60] 
(filled symbols) versus the mass accretion rate. The typical un- 
certainty for the colours is plotted close to the bottom right cor- 
ner. Triangles are upper limits for the accretion rates. 



have a disk fraction of 20%, while lower mass objects (2-0.8 
M ) have a disk fraction of 50%. 

3.2. Accretion rate and spectral energy distribution 

Figure [2] shows the distribution of mass accretion rates for the 
whole sample, as well as for the M01 group I and group II 
sources. The typical - median - accretion rate for each group 
is ~ 10~ 7 Mq yr -1 , the same as for the whole sample. Group 
II stars show accretion rates distributed roughly homogeneously 
around the median. Group I sources show more scatter, but 
their accretion rates also cover the full possible range for the 
HAeBe regime. Three group I sources, namely R Mon, VY Mon 
and LkHa 234, are hot -T* > 12000 K- HAeBes with strong 
Balmer excesses most probably associated with high accretion 
rates (Paper II). Even considering these objects, the null hypoth- 
esis that both M01 groups are drawn from the same accretion rate 
distribution is supported by a two-sample Kolmogorov-Smirnov 
(K-S) test with a probability of 75%. 

A larger number of M01 group II stars, compared to group 
I, is found in all works dealing with this classification (see Sect. 
12.1.11 and references below). The origin of this division could 
be understood by arguing that objects that show flared group I 
SEDs are in a younger, faster evolving stage than self-shadowed 
group II stars. This view, which would point to an evolution from 
group I to gro up II sources, is t heoretically and observation- 
ally supported (Ipullemondll2002t iDullemond & Dominikll2004t 
lAckeetalJ 12001 " It would imply that group I sources should 
show higher accretion rates than group II stars. As we men- 
tioned, such a conclusion cannot be drawn from our sample, and 
the implications deserve more study. On the other hand, our data 
agree with pre vious works that re ported larger dust grain sizes in 
group II stars (Acke et al. 2004]) - our median {3 value obtained 
from the fits in Fig. IA.2l is 0.6 for this group, against 1 .4 for group 
I sources -, and with works that reporte d that UXOrs tend to be 
group II sources (IDullemond et al.ll2003l) - LkHa 234 is the only 



UXOr (Oudmaijeretal ] l2001r.lRodgersll2003l) in our sample not 
belonging to that group, but this should be taken with caution 
given that its IR photometry is likely contaminated by the en- 
vi ronment (lEiroa et al.lll9 98) -. In addition to the interpretation 
in lDullemond et all d2003h . the fact that group II sourc es include 
most UXOrs could be related to an inclination effect (Natt a et alj 
119971: iNatta & Whitnevll2000|). which in turn could be affecting 
the SE P shape classification (Robitaille et al. 2007; Meiie r et alj 
I2008I 

From the observational perspective, the absence of a clear 
correlation between the accretion rates and the M01 groups was 
somehow expected. The M01 classification scheme strongly de- 
pends on the shape of the SED at long wavelengt hs (in particular 
on the IRAS 12, 25 and 60 fim fluxes; see van Boe kel et al . 2003, 
and Sect. 12.1. Ti l, and Fig.[3]shows that the mass accretion rate is 
not correlated with the colours at these bands. 

The contribution from the accretion shock is stronger at 
shorter wavelengths (see Paper II and referenc es therein), and 
could be measurable up to nIR wavelengths (Muz erolle et al.l 
2004). Figure |4] shows the intrinsic colour excess E(H - K)o 
against the mass accretion rates, along with the correspond- 
ing values for TT stars. Both parameters are correlated, with a 
Spearman probability of false correlation of 1.5 x 10~ 5 % for 
our sample. However, this relation becomes not statistically sig- 
nificant for M acc < 1.5 x 10~ 7 M0 yr -1 . The slope for stronger 
accretion rates is 0.12 + 0.05, e qual within the unc ertainties to 
that for CTT stars (0.10 + 0.05: lMever et al.lll997l) . Both sam- 
ples can be fitted simultaneously with a linear expression with 
slope 0.15 + 0.03. 

The similar correlation between the mass accretion rates and 
the H-K colour excesses followed by T Tauri and HAeBe stars 
points to a common origin for both trends, suggesting that the in- 
ner dust of HAeBe stars is not only heated by re-processing light 
from the central star, but also by the contribution of accretion. 
Because the relation between the mass accretion rate and the nIR 
colour breaks for M acc < 1.5 x 10~ 7 M0 yr -1 , this value can be 
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Fig. 4. Intrinsic colour excess in H - K versus the mass ac- 
cretion rate. The typical uncertainty for the colour is plotted 
close to the bottom right corner. Triangles are upper limits for 
the accretion rates. Crosses represent classical T-Tauri (bold- 
face crosses ) and weak T Tauri (normal crosses) stars from 
iMever et all (Il997l) . The best fits are E(H - K) = 1.51 (+0.26) 
+ 0.12(+0.05) x log M acc (dashed line, Pearson correlation co- 
efficient 0.50) for HAeBes with M acc > 1.5 x 10~ 7 M© yr _1 , 
E(H-K) = 1.26(±0.34) + 0.10(±0.05) x log M acc (dotted line, 
Pearson correlation coefficient 0.50) for CTTs, and E(H - K)o = 
1.63(±0.17) + 0.15(±0.03) x log M acc (solid line, Pearson corre- 
lation coefficient 0.65) for both HAeBes showing high accretion 
rates and CTTs. 



considered as an approximat e limit to divide passiv e and active 
disks in HAeBe stars (see also lvan d en Ancker 120051) . This is one 
order of magnitude larger than the corresponding accretion rate 
that divides weak and classical T Tauri stars (Fig. |4}, as can be 
estimated by assuming viscous accretion dissipatio n and stellar 
irradi ation as the only heating sources of the disk (ICalvet et al.l 
2000): 



M pa ~ 2 x 1Q-* M Q yr 



'( 



4000/:/ \2RqJ \0.5M o ) 



(2) 



where M pa , the accretion rate dividing passive and active disks, 
is derived assuming that the mass accretion rate through the in- 
ner disk is constant and equal to the rate at which mass is trans- 
ferred onto the central object. Applying Eq. |2]to a typical HAe 
star (T» = 9000 K, M» = 3 M©, R„ = 3 R©), the observational 
value 1.5 x 10~ 7 Mq yr _I is also recovered. This corresponds 
to log (L flcc /L Q ) ~ 0.60, which, along with the empirical rela- 
tions with the emission line luminosities provided in Paper II, 
indicates that active disks around HAeBes tend to show observa- 
tional thresholds log (L H JL Q ) > -1.5, log (L [0 ij63Oo/L ) > -3.7, 
and log (L firr /L ) > -3.2. 

Figure [5] shows the distribution of mass accretion rates ac- 
cording to the classification from the wavelength where the IR 
excess starts (groups JH and K). In this case, the K-S test pro- 
vides a probability of only 1.8% that both groups follow the 
same accretion rate distribution (null-hypothesis), against the al- 
ternative that both groups show different accretion rate distri- 




Fig. 5. Histogram showing the mass accretion rates for the 38 
stars in the sample (shaded bars), and the stars belonging to the 
groups JH and K (half of the stars in each group; violet and red 
lines). The probability that both groups are drawn from the same 
accretion rate distribution is < 0.018. 



butions. Moreover, if the alternative hypothesis is that the stars 
in the group JH have a smaller cumulative distribution function 
than the remaining objects, then the null hypothesis is rejected 
with a probability of 99.1 %. Therefore, we conclude that there 
is a significant trend for the stars with signs of inner dust dissi- 
pation to have lower accretion rates (< 10~ 7 Mq yr~' ), being the 
objects with IR excesses starting at / or H the strongest accre- 
tors. The typical - median - mass accretion rate for these stars is 
~ 10~ 6 Mq yr~\ one order of magnitude higher than that for the 
group K objects. This difference reaches two orders of magni- 
tude if mean values are adopted instead of medians. Some stars 
deviate from the general trend that relates the accretion rates 
and the wavelength where the IR excess starts. For instance, HD 
141569 shows an SED typical of a transitional disk (see the cor- 
responding panel in Fig. IA. Il l, pointing to significant inner dust 
dissipation and/or grain growth. At the same time, there are no 
signs of significant inner gas dissipation, since its mass accretion 
rate is only slightly lower than the median value for the sample 
(1.3 x 10~ 7 against 2.4 x 10~ 7 Mq yr 1 ). 

Finally, the median /3 value obtained from the fits of the mm 
fluxes (Fig. IA.2l ) has a similar value, about 1.1, for both groups 
JH and K, which contrasts with the results mentioned above 
about the differences for the M01 groups. This suggests that a 
lower nIR emission is is not necessarily linked to grain growth as 
determined by mm measurements; these sample the outer disk, 
and are more closely related to the possible decrease of the disk 
flaring angle associated to M01 group II sources. 

3.3. Accretion rate and disk mass 

Figure|6]shows the mass accretion rate against the disk mass for 
the stars in our sample, as well as for T Tauris. The Spearman 
probability of false correlation between log M acc and log 
is 2.9 % for our sample, which is best fitted from log M acc = 
1.1(±0.3) x log Mdisk - 5.0(±0.5), with a Pearson correlation 
coefficient of 0.60. Most stars in the sample deviate from that 
expression by less than one order of magnitude (dashed lines in 
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Fig. 6. Ma ss accretion rate against disk mass . Crosses are T Tauri 
stars from lHartmann et all (11998b and lNaiita etail (120071) . Dark 
crosses are transitional disks. The values for our 26 HAeBes with 
disk mass estimates are represented with solid symbols, triangles 
for the upper limits. Group JH stars are plotted in violet and 
group K stars in red. The solid line represents the best fit to our 
sample and the dashed lines + 1 dex. The close binaries for both 
T Tauri and HAeBe stars are surrounded by solid circles. 



Fig. [6]). The slope for T Tauris is also close to unit y, once they 
are divided into transition and non-transition disks (Naiit aet alJ 



2007), which is t he value expected from simple visco us dissipa- 
tion models (see Hartmann et al.|[l998l: iDullemond et alj|2006t 
Naiit a et alJ 12007 ). The trend could break for HBe stars more 
massive than those studied here, given the uncertainty of their 
accretion rates, and because some of their disk masse s are lower 
than those for T Tauri stars (lAlonso-Albi et al.ll2009l) . 

Given that M disk (r')/M disk (r) = (r'+l)/(r+l), with M disk (r) 
the total disk mass from Eq. Q~]for a gas-to-dust ratio r, a different 
value f would shift the data in the x-axis of Fig. [6] For instance, 
if f — 10, instead of the assumed value r = 100, then the x values 
would shift 1 dex to the left but the M acc -Mrf, S A relation would 
still have a slope close to 1 . In conclusion, the slope of this trend 
should not be affected by the chosen value of the gas-to-dust 
ratio, if this is roughly the same for all objects. 

Our sample is divided in Fig.[6]according to the wavelength 
where the IR excess starts, and considering the close binaries 
separated by less than 1000 AU, in analogy with the study in 
Naiita et aTT d2007|) for lower-mass stars in Taurus. In this case 
it was reported that stars with transition disks, when compared 
with non-transitional ones, tend to show accretion rates ten times 
lower for a given disk mass, and median disk masses about 
four times larger. In addition, close binary stars were reported 
to show slightly lower median disk masses than non-binaries. 
Although most of our group K stars are too far from being tran- 



sitionals, from the definition considering the wavelength where 
the IR exc ess is apparent, tw o of the three mentioned conclu- 
sions from lNaiita et al.l (120071) also apply to the HAeBes. First, 
the mass accretion rates tend to be lower by a factor > 10 for 
our stars with IR excesses starting at longer wavelengths (Sect. 
13. 2\ . Second, the close binaries in our sample have lower me- 
dian disk masses compared to systems without such compan- 
ions (1.2 x 10~ 2 a gainst 3.1 x 1 0~ 2 M ^). However, in contrast 
with the results in Naii ta et ail (j2007), the disk masses of the 
HAeBes with signs of inner disk dissipation are typically three 
times lower (median values) than those with IR excesses mea- 
sured at the J or H bands - considering the uncertainties in Table 
[1] the mean disk mass for group JH is 0.24 ± 0.02 M Q , while 
that for group K is only 0.066 ± 0.003 M . This difference be- 
tween TTs and HAeBes could be pointing to different physical 
mechanisms dominating inner disk clearing in both regimes (see 
the following section). 

4. Discussion 

The similarities found when the accretion rate is related to the 
nIR colours (Sect. I3.2l i. and disk masses (Sect. [33}, point to a 
common origin that explains these trends for both the T Tauri 
and the intermediate-mass regime. The relation with the nIR 
colour excess extends the corresponding one for TTs to HAeBe 
stars (Fig. 3]l. Both trends can be understood if the dust of the in- 
ner disk is heated not only by reproces sed stellar light, bu t also 
by the viscous accretion contribution (Cal vet et al.l l2000). The 
slope relating the accretion rate and the disk mass is close to 
unity (Fig. [SV which is also expec ted from simple viscous disk 
models (e.g iDullemond et alj|2006l) . 

Apart from the analogies mentioned, we found two major 
differences when our data are related to the corresponding for T 
Tauri stars. First, our results are consistent with a slightly shorter 
inner gas dissipation timescale (Sect. [XT). As argued in that sec- 
tion, this could have implications on the physical mechanism 
able to form planets around stars more massive than the sun on 
such a shorter timescale. Second, the relative position of HAeBe 
stars with signs of inner dust dissipation in the ld a0 c-^disk plane 
(Fig. |6j, compared to the remaining HAeBes without signs of 
disk clearing, differs from a similar comparison between transi- 
tion and classical disks around TTs. In particular, our group K 
sources show lower accretion rates and disk masses than group 
JH sources, whereas transitional T Tauri stars show lower accre- 
tion r ates but larger disk masses than classical TTs (Naiit a et al.l 
2007). The M acc -Md,-^ plane is a potential tool for studying the 
physical mechanisms that dominate circumstellar disk dissipa- 
tion. In summary (see lAlexander & Arm itage 2007; Naiit a et al.l 
2007, for more details), large disk masses and low mass accre- 
tion rates, compared to those typically shown by non-transitional 
disks, are expected if Jovian planets are formed around stars with 
signs of inner dust dissipation. Large disk masses are necessary 
to form this type of planets, and at the same time, the inner gas 
content - and therefore the accretion rate - is strongly dimin- 
ished since pail of the gas otherwise destined to the central star 
accretes onto the planets. In contrast, photoevaporated disks are 
in principle more consistent with both low disk masses and ac- 
cretion rates. It is unknown if this type of c onclusions from the 
modelling in lAlexander & Armitagel (12007b could be applied to 
stellar masses higher than 1M Q . As we mentioned (see also Sect. 
13.31 ). most of our stars with signs of inner dust dissipation have 
mass accretion rates and disk masses lower than the correspond- 
ing median values for the remaining sources by a typical fac- 
tor of ~ 10 and 3, respectively. Some examples are HD 34282, 
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HD 34700, HD 142666 or BF Ori, as well as transition disks 
associated with distant and close binary stars (HD 34700 and 
HD 141569, respectively). The lower values for both M acc and 
M-disk could be expected from photoevapor ation sweeping out 
both the circumstellar g as and dust content (IClarke et al.l l2001; 
lAlonso-Albi et al.1 12009). Indeed, from the theoretical point of 
view photoevaporation could be more relevant for dissipating 
disks in the HAeBe regime than in the T Tauri one (see e.g. 
iTakeuchi et al.ll2005l) . On the other hand, T Ori is the only star 
in our sample with signs of inner dust dissipation that shows a 
lower mass accretion rate and a higher disk mass than the corre- 
sponding median values for t he stars with an IR excess starting 
at J or H bands. As argued bv lNaiita et al.l (120071) for the transi- 
tion disks in Taurus, those properties are anticipated by several 
planet formation theories, suggesting that Jovian planets could 
be playing a significant role in explaining the inner disk dissi- 
pation of T Ori. The possible presence of a T- Tauri type stellar 
companion located at ~ 3500 AU, considering a distance to th e 
system of 470 pc (iLeinertet alJll997t iMontesinos et al.l l2009h . 
could also be affecting the position of T Ori in the M acc -M c usk 
plane. CQ Tau shows similar characteristics to T Ori, with signs 
of inner disk dissipation, a low value of M acc and a high M^, 
compared with the medians for the whole sample. Thus, the for- 
mation of Jovian planets could also be playing a role for this star. 
Grain growth is also apparent from the fitted value of /3 (~ 0.6). 

In conclusion, we have provided several observational in- 
dications that could be consistent with various disk dissipation 
mechanisms - viscous accretion, photoevaporation, planet for- 
mation, grain growth -. The dominating mechanism could differ 
depending on each object. To establish whether there could be 
a dominating disk dissipation mechanism for the whole HAeBe 
regime would require additional observational and theoretical ef- 
fort. 

The discussion above must take into account that disk 
masses estimated from mm measurements could be affected 
by strong uncertain t ies. F or a sample of T Tauri stars, 
lAndrews & W illiams ( 20071) found that disk masses from mm 
continuum emission underestimate those fro m accretion by 
roughl y one order of magnitude. Following lHartmann et al.l 
(1998), the amount of material accreted from to to t^s - with 
to the age of the star and t^s the time to reach the main sequence 
- is a lower limit for the disk mass at to, and is given by 



M, 



disk 



f 



M acc (t)dt. 



(3) 



Assuming M acc (t) = AV n + B, with rj > 1, A, B constants that can 
be determined from a measurement of the mass accretion rate at 
a given to, and under the condition M dcc (tMs) ~ 0, we derive 



M disk > 



( M. dCC (t ) ) 



W - 1, 



x 



1 
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(5) 



v l l MS ' 

which for t^s » to, yields 

^ M. dcc (t ) 

M disk > — x f . 

rj- 1 

Equation |5] was applie d by lHartmann et al.l (1 19981) and 
Andrews & Williams (2007) to T- Tauri stars. It gives practically 
equal values to Eq. |4] for these objects, and for many of our 
stars, but Eq. [4] provides disk masses that can be one order of 
magnitude larger than those from Eq. [5] for the most massive 
stars in our sample. These objects are expected to evolve very 




log M .. 



-1 1 

| (from mm measurements) 



Fig. 7. Disk masses from accretion versus those from dust emis- 
sion at mm waveleng t hs. Cr osses are T Tauri stars taken from 
lAndrews & Williams! (|2007 ) and circles the estimates for our 
sample. The solid line represents equal values and the dotted line 
bisects both samples, indicating that disk masses from accretion 
are typically one order of magnitude larger than those from mm 
measurements. The filled circles are derived from Eq. |5]and ^ 
= 1.8. The error bar at the right side represents the uncertainty 
coming from using 77 = 1.1 (upper limit) or 3.2 (lower limit). The 
open-dotted circles are the estimates from Eq. |4] 7 7 _ 1.8, and th e 
theoretical times to reach the main sequence 



}.|4J 77 = 1.8, anatne 
from TavTerl(ri99l . 



fast because their ages are similar to the theoretical times to 
reach th e main seque nce (i.e. to ~ Ws, see Fig. 4 in paper II 
and e.g. lTavlerlll994 . Figure [7] compares the disk masses de- 
rived from the mass accretion rate to those from the cold dust 
emission. The results for our sample extend the findings for T 
Tauri stars, indicating that disk masses inferred from mm mea- 
surements are smaller than those from accretion by an order o f 
magnitude on average. Following lAndrews & Williamsl (120071) . 
this may be the result of particle growth (up to approximately 
meter sizes), which leads to an overestimate of the opacity used 
to compute from the millimeter luminosity. This view is 
supported by the low /3 values derived for most stars in our sam- 
ple with SED fitting at mm and sub-mm wavelengths (Fig. | A.2t . 
which are similar to those for T Tauris dRicci et al.ll2010l) and 
point to substantial grain growth in the CS environm ent of pre- 
main sequence stars (see lAndrews & W illiams 2007, and refer- 
ences therein). For a fixed value of the dust opacity, the use of 
Eq.Q]with a gas-to-dust ratio ~ 1000 yields disk mass estimates 
an order of magnitude higher, which would roughly match those 
derived from accretion. Although there is no strong evidence for 
such a high ratio compared to what is typical for the interstellar 
medi um, this possibility cannot be excluded in som e objects (see 
e.g.lMartin-Za idi et al.ll2007T: ICarmona et al.ll2008l) . 
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5. Summary and conclusions 

We have related the mass accretion rates with the stellar ages, 
SED properties, and disk masses of a sample of 38 F-G, HAe 
and late-type HBe PMS objects. We have compared the trends 
with the corresponding for lower mass T Tauri stars. Our main 
results and conclusions are: 

- The mass accretion rate decreases with the stellar age. This 
decline can be fitted to M acc (t) oc exp(-f/r) or to M acc (t) oc 
l~ n . The first expression provides a dissipation timescale t = 

Myr. The second expression is the same as used in 
simple viscous dissipation models for CTT stars, and yields 
j] = 1.8+07. T ne lar 8 e error bars reflect intrinsic scatter, un- 
certainties and upper-lower limit estimates for the accretion 
rates and stellar ages. The values obtained should be viewed 
with caution, given that both the the stellar age and the ac- 
cretion rate also depend on the stellar mass. However, our 
expressions point to a slightly faster dissipation of the inner 
gas in HAeBes, when compared to lower mass T Tauri stars. 

- Our sample does not show any apparent correlation between 
M acc and the SED shape according to the Meeus et al. (2001) 
classification scheme. This is more closely related to the 
grain size or the UXOr behaviour, as was previously re- 
ported. 

- The correlation between the near-IR excess and the accretion 
rate extends from the CTT to the intermediate-mass regime. 
The mass accretion rate that divides active and passive disks 
rises one order of magnitude for our sample (~ 2 x 10~ 7 
M0 yr~'), as can be inferred assuming viscous accretion and 
stellar irradiation as the only heating sources of the inner 
disk. 

- In agreement with simple viscous disk models, the mass 
accretion rate correlates with the disk mass as M acc oc 
^idisk ±03 \ with a scatter lower than + 1 dex for most 
of the stars, which is roughly followed by CTT stars once 
they are classified into transition and non-transition disks 
dNaiita et al.ll2007l) . 

- Most stars in our sample with signs of inner dust disk clear- 
ing (IR excesses starting at the K band or longer wave- 
lengths) show lower accretion rates and disk masses than 
those for the remaining objects. This contrasts with previ- 
ous results for transitional T Tauri stars, and could suggest 
that a different physical process - such as photoevaporation - 
plays a major role in dissipating disks in HAeBes. 

Finally, we extended the previously reported mismatch be- 
tween the disk masses derived from cold dust emission and ac- 
cretion to the HAeBe regime. More efforts are necessary to break 
that inconsistency, which is most probably related to the way that 
disk masses are estimated from mm measurements. 
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Appendix A: Tables with photometry and SEDs 



Table A.l. Ultraviolet, optical and near-IR photometry 
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NotesandreferencestoTablelA.il Numbers are in magnitudes. Most UBVRIJHK photometry is taken from the 
simultaneous EXPORT multi-epoch observations ( Oudm aiier et al] 1200 ll : lEiroa et al ] 120021) . selecting the data at 
the brightest V-magnitude (italic n umbers for the ma g nitudes not measure d simultaneously to V) . *SIMBAD 
(http://simbad.u-strasbg.fr/simbad/), 'Malf ait et all dl998h . ^Coulson et al.Ul998h. ^Sylvester et all dl996l) (inc l udes L' 
magni tu des forHD 141569, H P 142666 and HP 14443 2) jHillenbrand et al]dl992l) jBigav & Garniejj 19701), iMerfnl 
(2004), iHerbst et al.l d!982h . iGlass & PenstorJ dl974h . ' iKenvon & Hartmannl dl995h . " iMonnier et alj d2009h . Good 
quality IUE spectra are available for the stars indicated in Col. 2, which were used to derive synthetic UV photometry 
to fit the SEPs. Uncertainties are available in the references, and are typically < 0.05 magnitudes for the EXPORT data. 
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Table A.2. Mid-infrared photometry 



Star 


ISO 


F, 


Fl2 


Fl8 


F 25 


F60 


F 65 


F90 


F100 


Fl4o(Fl6o) 


Add. Refs 


HD 31648 




9.118 


12.25 


8.070 


10.34 


11.11 


10.17 


11.73 


12.52 


11.85 


(1) 


HD 34282 




0.736 


0.700 




1.630 


10.40 


10.04 




10.72 


9.745(9.353) 
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V 




0.600 
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<7.360 






HD 141569 


V 
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V 
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V 
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HD150193 


V 


13.00 
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5.839 
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< 16.25 




(1,2) 
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V 


15.02 


18.20 




20.99 


28.24 






<40.62 




(1,2,3) 


HD 179218 


V 


14.74 


23.44 


27.88 


43.63 


29.92 


19.43 


16.80 


17.35 


9.263 


(1) 


HD 190073 




5.734 


7.160 


4.457 


5.530 


1.920 




1.119 


< 1.000 




(1) 


AS 442 






<5.188 




2.690 


<9.620 






<7.215 






VXCas 




0.7554 


1.540 


1.327 


2.640 


1.400 




0.7453 


<15.17 
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0.481 


0.475 


0.767 


0.116 


0.138 
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<0.257 






BOCep 




0.084 


<0.070 


0.118 


0.285 


1.430 




1.566 


<3.770 
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1.858 
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V 


38.73 


72.80 


79.53 


119.0 


474.0 






880.0 




(2,4) 
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34.07 


54.65 


79.24 


132.1 


121.3 


129.5 


78.13 


148.5 


66.37(62.26) 


(1,2) 


VY Mon 






42.20 


56.22 


78.50 


133.0 


129.0 


115.1 


260.0 


100.1 


(4,5) 


51 Oph 


V 




15.67 


9.241 


10.19 


1.060 






<5.970 




(1) 


KK Oph 




9.100 


9.870 


8.276 


9.560 


6.140 


3.905 


4.431 


<17.75 




(2) 


TOri 






















(2) 


BFOri 


V 




0.960 




0.800 


<2.460 






<29.81 




(2) 


CO Ori 




1.677 


1.580 


1.537 


1.400 


<2.030 






<10.88 






HK Ori 




2.708 


3.800 


2.882 


4.080 


< 1.640 






<70.37 




(2) 


NVOri 




1.247 




















RYOri 




0.665 


0.750 


0.788 


0.710 


<2.280 






< 19.09 






UX Ori 


V 




2.680 


3.112 


3.690 


2.850 


2.736 


2.129 


3.760 


4.282 


(2,6) 


V346 Ori 






0.310 


0.590 


1.250 


5.630 


4.054 


3.793 


3.470 


3.917 


(7) 


V350 Ori 




0.556 


0.940 


1.008 


1.750 


1.190 






<21.94 






XY Per 




4.302 


3.850 


3.404 


4.060 


4.910 


3.958 


4.367 


< 10.27 


5.323 




VV Ser 


V 


4.261 


4.610 


3.259 


3.770 


6.340 






24.29 




(2,8) 


CQTau 


V 


4.855 


6.210 


12.88 


20.66 


21.88 


15.89 


16.25 


13.65 


10.30 


(6) 


RRTau 


V 


1.365 


1.740 


1.367 


2.200 


4.460 






36.68 




(2) 


RYTau 




12.28 


17.47 


15.43 


26.07 


15.29 






13.65 




(2,6,9) 


PX Vul 




0.433 


0.550 


0.569 


0.860 


<0.460 






<5.060 






WWVul 


V 


1.502 


1.660 


1.786 


2.210 


1.560 




1.357 


<9.590 




(6) 


LkHa 234 


V 


5.844 


14.78 


18.72 


78.96 


687.7 






1215 


757.2(880.6) 


(1,4) 



Notes and references to TablelA.2l Fluxes are m Jy. IRAS fluxes at 12^ 25^ 60 and 100 /7m are from S 1MB AD 
(http://simbad.u-strasbg.fr/simbad/) by default, and from the IRAS Faint S ource Reject Catalog (I PAC 1992) for AS 442, IRAS 
Faint Source Catalog, |b| > 10, Version 2.0 for BH Cep and BO Cep, and from lCasev & Harperld 19901) for VY Mon. AKARI fluxes at 
9, 18, 65, 90, 140 and 160 /im are from AKARI/IRC All-Sky Survey Point Source Catalogue (Version 1.0) and AKARI/FIS All-Sky 
Survey Bright Source Cat alogue (Version 1.0). Additional references are given in Col. 12 for photometry in the following wave- 
lengths: (1): 10.7 um from lMonnier et all d2009l) : (2): 10.5(N) and 20. 1(0) um from lffllienbrand et all dl 992): (3): Spitzer IRAC 5.8 
um from the GLI MPSE catalogue; (4): Spitzer IRAC 3.6, 4.5, 5.8 and 8 . um from the "Spitzer survey of young stellar clusters" 
dGutermufh e7aTl 2009): (5): 10.5(N) and 20.1(Q) fim from lHerbst et all d!982l): (6): ISO photo metry at 6.9, 9.6, 17, and 28.2 fim 
from lThi et al.l d200lh : (7): Spitzer IRAC 5.8, 8.0 and MIPs 24 um from lHernandez et al] O0 06); (8,9): Spi tzer IRAC 3.6, 4 5, 5.8 , 
8.0 and MIPs 24, 70 fim from the "Spitzer survey of Serpens YS O population " dHarvev* etail2 007) and from Robit aille et all d2007l) . 
respectively. Photometry extracted from ISO data is available in Merin (2004) for the stars indicated in Col. 2. The typical uncertainty 
is < 10 % for IRAS and <, 5 % for AKARI data. 
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Table A.3. Sub-mm and mm photometry for a sub-sample of 
stars 



Star 


F( 1.3mm) 


F(A) 


Ref. 




(mJy) 


(mjy)((mm)) 




HD 31648 


360±24 


3300(0.45); 780(0.85); 235(1.4); 39.9(2.7); 35.2(2.8) 


(1,2,3) 


HD 34282 


110+10 


2700(0.45); 409(0.80); 360(0.85); 183(1.1); 24(2.6); 5.0(3.4) 


(4,5,6) 


HD 34700 


11.7+1.1* 


218(0.45); 41(0.85); 39(1.1) 


(6,7,8) 


HD 141569 


9.0±4.0 


65(0.45); 11(0.85) 


(5,9) 


HD 142666 


127+9 


1090(0.45); 351(0.80); 180(1.1); <63(2.0) 


(6) 


HD 144432 


37±3 


103(0.80); 129(0.85); 69(1.1) 


(5.6.9) 


HD150193 


45±4.5 


101(0.85) 


(5,10) 


HD 163296 


780+31 


10660(0.35); 5560(0.45); 3380(0.60); 2650(0.75); 2320(0.80); 1920(0.85); 1020(1.1) 


(11) 


HD 179218 


71±7 


7.6(2.6) 


(12) 


VXCas 


<6 


... 


(13) 


SV Cep 


7±2 




(13) 


V1686Cyg 


<4.8* 




(14) 


RMon 


100+10 


1980(0.35); 1560(0.45); 217(0.80); 77(1.1); 1.17(20); 0.37(36); 0.29(60) 


(11,15,16) 


VY Mon 


120+12 


18000(0.37); 1.5(60) 


(10,17) 


KKOph 


52+5.2 


109(0.80); 36(1.1); 0.24(36) 


(5,16) 


TOri 


130+28 




(15) 


BF On 


6+2 


... 


(13) 


CO Ori 


<4.0 




(13) 


HK Ori 


<100 




(15) 


UXOri 


23+2 


3.8(2.6); 0.8(7.0) 


(13,18,19) 


VVSer 


<20 




(15) 


CQTau 


143+8.4 


2.6(7.0) 


(12,19) 


RRTau 


<20 




(15) 


RYTau 


212+21* 


2439(0.35); 1920(0.45); 560(0.85) 


(20,21) 


WWVul 


10.5+1.0 




(13) 


LkHa 234 


827+82.7* 


60000(0.45); 4200(0.80); 3400(0.85); 2.67(20); 1.82(36); 1.65(60) 


(16) 



Notes. The second column shows the fluxes used to compute the disk masses in Table [T] Those tagged with * were measured at 
1.35 mm (HD 34700), 2.70 mm (V168 6 Cyg) and 1.20 mm (RY T au an d LkHa 234). An uncertaint y of 10% is assigned w hen 
that is not provide d in t he reference. (l):lMannings & Sargentl dl997l). (2 ): ISandell & Weintraubl (120031) . ( 3): iPietu et al.1 d2006l) .(4): 



Pietu et all d2003l) (5): [Sandell et all d201 ll). (6): ISvlyester et al.l d 19961) (7)TISvlvester et al.l J200lh. (8): ISheret et all d2004h. (9): 
Walker & Butneri(ll995h . (10):|Alonso-Albi et al.l(T2009h (1 l):lManningsl d!994h. (12):lMannings & Sargentl d200rJ).lT3 ) :lNatta et al.1 
dl997l). (14):lNatta et alJ d200o|) . (15): lHillenbrand et al Jdl992h7(16):lPezzuto & StrafellaUl997t). (17):ICasev & Harpe3 ( ll990h . (18): 
lAcke & van den Anckej mM) , (19): lTesti et al.l d2001h . (20): lRobitaille et all j2007Tr (21)-TAltenhoff et al.| d 19941) 
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Fig. A.l. SEDs of the stars in the sample. Triangles are photometric upper limits. For each panel, the best fit is the black-solid line 
and th e models fitting with^ 2 -)A est ^ 3n are the grey-solid lines, where n is the number of datapoints per star (see Rob itaille et al.l 
I2007I) . The best stellar (de-redenned) photosphere is indicated with a dashed line. 
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Fig. A.2. Subsample of stars with fluxes available at A > 0.35 mm (log v [GHz] < 2.9, see Table lA.31 . The best graybody fits yielding 
the indicated values for To and (3 are shown. Fluxes at A < 2.0 mm (log v [GHz] > 2.1) are weighted stronger in the fits. 



21 



